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In the framework of Horava-Lifshits theory, we study the eigenvalues associated to the Wheeler- 
DeWitt equation with the meaning of vacuum states, i.e. cosmological constants. The explicit 
calculation is performed with the help of a variational procedure with trial wave functionals of the 
Gaussian type. We analyze both the case of the detailed balanced condition and the case without 
it. We find the existence of an eigenvalue depending on the coupling constant and on the physical 
scale 
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I. INTRODUCTION 

o 

■ The cosmological constant represents one of the challenges of this centuryfT'|. If one is tempted to compute it from 
the zero point energy of some physical field, one discovers an enormous discrepancy with the observation data. This 
discrepancy amounts to be of 10^^° order of magnitude: this is known as the cosmological constant problem. Many 
^—y' attempts to explain such a discrepancy have been done, but they appear to be far from being satisfactory. Recently, 
O ; Hofava proposed a modification of Einstein gravity motivated by the Lifshitz theory in solid state physics [2] [3]. 
^ '' Such modification allows the theory to be power-counting ultraviolet (UV)-renormalizable and should recover general 
$H relativity in the infrared (IR) limit. Nevertheless Hofava-Lifshitz (HL) theory is non-covariant. Indeed, in this 
approach space and time exhibit Lifshitz scale invariance of the form 
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t tt and ^ fe* (1) 



with z > 1. z is called the dynamical critical exponent and in the present case it is fixed to z = 3. The breaking of the 
I 4D diffeomorphism invariance allows a different treatment of the kinetic and potential terms for the metric: from one 
side the kinetic term is quadratic in time derivatives of the metric, form the other side the potential has high-order 
space derivatives. In particular the UV behavior is dominated by the square of the Cotton tensor of the 3Z? geometry 
, by means of a contribution to the propagator leading to a renormalizable power-counting theory. The original HL 
' theory is based on two assumptions - detailed balance and projectabilityQ^. The projectability condition is a weak 
^\ \ version of the invariance with respect to time reparametrization and therefore to the Wheeler-DeWitt equation[T7|. 
■ Motivated by these interesting features, we ask ourselves if the problem of the cosmological constant in HL theory 
has better chances to be solved than in ordinary Einstein gravity^. Indeed, the modification of the gravitational field 
• ^ at short distances could produce new contributions that allow what in ordinary Einstein gravity is forbidden. In 
ordinary Einstein gravity, there exists a well accepted model connected with the cosmological constant: this is the 
^ Friedmann-Robertson- Walker metric, whose line element is 

ds^ = -N'^dt^ + (t) dnj. (2) 

dr^l is the usual line element on the three sphere and is the lapse function. In this background, we have simply 
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For simplicity we consider the case in which k ~ 1. The generalization to fc = 0, — 1 is straightforward. The Wheeler- 
DeWitt equation for such a metric is 



i/^- (a) = 



92 9^2 ^ 2 A 4 



^ (a) = 0. (4) 
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1 Different aspects of HL theory are discussed in Refs. H 0, S II 13, [H [U, [12, El, Q El, El ■ 
^ See also Ref. llSll . for a different approach on the cosmological constant problem in HL theory. 



2 



It represents the quantum version of the invariance with respect to time reparametrization. If we define the following 
reference length ag = \/ 4 j the Eq. ^ assumes the familiar form of a one-dimensional Schrodinger equation for a 



particle moving in the potential 



Uia) = 



4G2 



(5) 



with total zero energy and without a factor ordering. Eq.(I3]) can be cast into the following form 



JVa^* (a) 
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J Pa** (a) 




4G2 



(6) 

which appears to be an expectation value. In particular, it can be interpreted as an eigenvalue equation with a weight 
factor on the normalization. The application of a variational procedure in ordinary gravity with a trial wave functional 
of the form 



* = exp (-/3a^) 



(7) 



shows that there is no real solution compatible with the procedure. The purpose of this paper is to obtain an eigenvalue 
equation for the cosmological constant like the one in Eq.B but in HL theory with the help of the Wheeler-DeWitt 
equation. Nevertheless, as pointed out by MukohyamaflSJ], there are four versions of the theory: with/without the 
detailed balance condition; and with/without the projectability condition. In this paper, we will consider the problem 
with and without detailed balanced condition. 

II. HL GRAVITY WITH DETAILED BALANCED CONDITION 

With the assumption of the detailed balanced condition, the action for the Hofava theory assumes the form 

S = I dtdPx{CK- Cp), (8) 

where 



(9) 



is the Lagrangian kinetic term. The extrinsic curvature is defined as 

1 



(10) 



with K = K^^ Qij is its trace and Ni is the shift function. For the FRW metric, there is no shift function. The 
Lagrangian potential term Cp is 



kV /1-4A 



where 



1 



is the Cotton tensor. Plugging the FRW metric ([2]) into Cp, we obtain 

3 6Aw 



Cp = 



3A2 



(l-3A)a4 (l-3A)a2 (1 - 3A) 



and the action for the potential part reduces to 
Sp^ 



dtd^xCp = ^ / dtN2TT^a^ 



Ex/ 



6A 



w 



3A2^ 



(l-3A)a4 (l-3A)a2 (1 - 3A) 



(11) 
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(14) 
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Concerning the kinetic part one gets 



CK=N^^{K'^K,,-XK^)=a^sm^Xsin0^[-] (1 - 3A) (15) 
and the corresponding action is 



Sk^ f dt(fxCK = I dt2n^a^— ( (1 - 3A) 
Jt,xi Ji N \aj 



(16) 



Now, we can compute the canonical momentum and we find 



TTa = = — (1-3A), (17) 

da K 

where we have set N = I. The rcsuhing Hamiltonian is computed by means of the usual Legendre transformation 
leading to 



H = (fxH = / (fx [iTaa - C] 



127r2a(l - 3A) 



a(l-3A) (1-3A) (1-3A) 



The classical constraint can be read off quite straightforwardly. We obtain 

ttI + 9^ V [l - 2Kwa^ + ^w^^] = 0. 



(18) 



(19) 



This is the classical constraint, which has been analyzed in Ref.[lO|. However, we are interested in the Wheeler-DeWitt 
equation associated to Eg. fig]) . This is simply given by 



[1 - 2Kwa^ + A^a"*] ^' = 0, 



(20) 



where we have promoted the canonical momentum tt^ to an operator tt^ = —ida- Even if the WDW equation can be 
solved. Aw cannot be treated as an eigenvalue, because Ayy is ent ang led with the scale factor. Things do not change 
if we make an analytic continuation of parameter fi i/x[l,[ld[il|,[l2|,[ll,[l3,[23|. To this purpose we need to slightly 



modify the action satisfying the detailed balanced condition j4!| by adding an appropriate IR term of the form 

^2 



L,p 



8(1 - 3A) 

where A is a dimensionless parameter. Then Ea. (|T8l) modifies into 

3 ^ 6{n'^~Aw)a 3A2j,a3 



H = 



2A/a3 
(1 - 3A) (1-3A) 



a(l - 3A) 



127r2a(l-3A) 4 
Simplifying somewhat, we obtain for the WDW equation 



(1 - 3A) 



da? 



Stt^H^ [3 + 6 (n^ ~ Aw) + 3A^a'* - 2Ap;'a^\ 0, 



which looks like Eq.ljl]). Eg. f^H]) can be recast into the form 



\Va^* (a) 



37rV' (3 + 6 iy? - Aw) a? + ZA\^a^) * (a) 



^T^^A^i' 



J Va^* (a) [a4] * (a) 

We compute the value of A, by adopting a variational technique with a trial wave functional of the form 

^ = exp (-/Ja^) . 



(21) 



(22) 



(23) 



(24) 



(25) 
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We find 



V 8/3 

which represents the normahzation of the wave function. The other terms become 



da* 



da** [2/? -4^V]* = 2/3^^-4/3^^1(1 



IStt^A^^ (m^ - Aiv) / da**a2* = 18n*^'^ (^^ _ 



TT /I 



8 V4 



and 



da**a''* = 



8 V16 



Phigging Eqs.lim . . {2M into Eq.dH]), one gets 



Simphfying somewhat, we obtain 



8 (l6^ 



8 \ le/^ 



27 



^TT^A^s = /33 + OttV'/?' + ^vrV' (/i' - Am.) /3 + ^A^^^^^^ 



We now compute the extreme of A = A (/3) to obtain 



whose solutions are 



^ = 
dp 



fii.2 = 



3/3^ + 187rV^/3 + (m^ - A^) = 0, 



SttV 



A 



We can recognize three cases: a) x = 1, b) a; 1 and c) < a: ^ 1. 

a) X — 1. This imphes that 

/3i,2 - -SttV (1 ± 1) ■ 
Both solutions lead to a non-normalizable wave-function and therefore will be discarded. 

b) a; ^ 1. For this choice, one solution must be discarded for normalization conditions and we remain with 



/3i ~ -3^V 1 - \/i + TT-r - ^ ^^ \hr 



67r^ 



(26) 

(27) 
(28) 

(29) 

(30) 

(31) 
(32) 

(33) 
(34) 
(35) 



Finally, we have the last case 
c) < a; <C 1 . The two solutions can be approximated with 



/3i,2 ^ -SttV ± 3M^7r%/ 1 



2„4 



67r4 



37rV (-1± 0.999) 



(36) 



and this is close to the case a) and therefore will be discarded. 



Therefore, only case b) can give potential solutions to the problem. Plugging solution ([35]) into Eq. (l3T|) . one 
approximately obtains 



A 



(37) 



Here the dominant contribution has the correct sign. Therefore we conclude that in this range there exists a family 
of solutions to the cosmological constant depending on the variable x. In the next section, we abandon the detailed 
balanced condition to see if we can obtain other different solutions. 
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III. HL GRAVITY WITHOUT DETAILED BALANCED CONDITION 



In Refs. [2l|, , a more general potential form which avoids the detailed balanced condition has been proposed. Its 
expression is 



+g5R {R''R.j) + 96R]RlR- + giRV^R + g^V ,RjkV' R'"" 



(38) 



Couplings (7a (a = . . . 8) are all dimcnsionless and powers of C, are necessary to maintain such a property of ga- 
Plugging the FRW background into £p, one gets 



9oC^ + 5iC^^^2^ + (6.92 + 33) + ^^6^ (36.g4 + 655 + ge) 



(39) 



The term goC^ plays the role of a cosmological constant. In order to make contact with the ordinary Einstein-Hilbert 
action in 3 + 1 dimensions, we set without loss of generality go(^ — 2 A and 51 = —1. In case one desires to study the 
negative cosmological constant, the identification trivially will be goC^ = — 2A . After having set N = 1, the Legendre 
transformation leads to 



Ti, = TTafl — Ck + Cp 



and the Hamiltonian becomes 

H = / (f xH = - 

where 



127r2a(3A- 1; 



+ 27r^a^ (i) 



2A 



6c 



a?{t) a^{t) a^{t) 



6ff2 + 9z = b 
36^4 + 655 + 56 = c 



The WDW equation can be easily extracted to give 



-2A 



6C^ 6C^6 6c 



^ = Q 



and adopting the same procedure of Eq. (j23p , we can write 



da^* lulm ^)o.„4„4, 



-247r*a* [t) 



-2A 



6c 



a2 (t) a4 (t) a6 {t) 



* = 0. 



(40) 



(41) 



(42) 



(43) 



(44) 



We can rearrange the previous expression to obtain 



da^* {nl^ + (^1447r4 [C^a^ (t) - C'b - ca'^ (t)] } * _ (3A - 1, 



487r*A. 



(45) 



It is now clear the role of A. As a trial wave functional, we adopt the same form of Eq. ([25|l and using the results of 
Eas. (P51 . . I 29p we can write 



/3V^+^144^' 











(^48.4A, 



(46) 



where we have used the following relation to compute 



c I da^ a = -c / — = 



exp(-2/3x) =-^r (-^^^c^j2^. 



(47) 
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By simplifying somewhat, one gets 



(1 + Ace) P^ + c 



-f3-P^C'^b 



where we have defined 



(3A-1) , 



By applying the variational procedure, one obtains 



dA 



3(1 + 45c) 



— = = 3(1 + 4Sc) (3^ - 25bCl3 + V 
dl3 4 

In order to have two real solutions, we impose that 



/?1,2 = 



> (1 + 4cc) -. 



(48) 



(49) 



(50) 



(51) 



Since the coefficients b and c are a linear combination of the coupling constants gt, for simplicity, we investigate two 

|; b) the case when = (1 + 4cc) |. 

2£bC 



cases: a) the case, when (1 + 4cc) |; b) the case when = (1 + 4cc) |. For the case a), we get 



Pi = 



3(1 + 45c) ■ 



Plugging Pi into Eq. ([i5|) . we obtain 



6(1 + 45c) 
If we approximate 1 + 45c ~ 45c, we get 



856^ 



9 (1 + 45c) 



= 9.^ i^A. 



^ 3c 



1 



26^ 
'9c 



= A<0, 



(52) 



(53) 



(54) 



because of the approximation chosen. Therefore this range of solutions will be discarded. For the case b) , we obtain 

/3i,2=/3=£. (55) 

Eq. (|48ll becomes 



126 



C«1 = A = ^. 
^ 3b 2 



(56) 



This implies 



36-^"- 



(57) 



IV. CONCLUSIONS 



In this letter, we have considered the recent proposal of Horava to compute the cosmological constant in a modified 
non-covariant gravity theory at the Lifshitz point 2 = 3. We have used the Wheeler-DeWitt equation as a backbone 
of all calculations with the cosmological constant regarded as an eigenvalue [23]. We have analyzed the situation with 
detailed balanced condition and without the detailed balanced condition. A variational approach with Gaussian wave 
functions has been used to do calculation in practice. The background geometry is a Friedmann-Robertson- Walker 
metric with k = 1. Concerning the detailed balanced condition, we have found an approximate eigenvalue not in the 
original configuration but with a slight modification in the IR region. The main reason is related to the fact that the 
original parameter A\y cannot be considered as an eigenvalue in this formulation. What we have found is that we have 
only results in agreement with a big cosmological constant, therefore with Planck era estimations. This means that 
the cosmological constant problem in this approach is not solved. It is interesting to note that only the dimensionless 
ratio Avy/z^^ comes into play. As regards the case without detailed balanced condition, we obtain two cases: 
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a) go is large. For example, it could be fine tuned to Planck era values. This implies that b = 652 +53 ^ 1- This 

can be achieved if 32153 ^ 1 or 6(72 — ~ff3- When g2,g3 <SC 1 we fall into a perturbative regime, while when 
652 — —93 this could not be the case. In this respect, this version of HL theory and the version with detailed 
balanced condition behave in the same way, except eventually for the smallness of the coupling constants 172 and 

53- 

b) go is small and fine tuned to go ^ 10^^^° as from observation. This implies that b — 632 + 53 S> 1, that it means 

that at least either g2 or g^ must be very large preventing a perturbative approach. 

Note that in the version of HL theory without balanced condition, there is no trace of the other coupling constants: 
only (72 and gs come into play. We can further observe that another big difference between HL theory with and without 
balanced condition is that the latter admits a possible transition from big values to small values of the cosmological 
constant. However, it is likely that an improvement could be obtained by introducing a renormalization group equation 
for (72 and (73 and a running A. Another improvement for having bounds on parameters can be extracted from the 
equation of state p ^ top with uj — —1. Using the Gaussian wave function's value of the exponent (5 obtained by the 
variational procedure into the equation of state, we can extract other informations on 52 and (73. Needless to say that 
the analysis has been done without any matter field and also this addition could improve the final result. 
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